Quantifying the impact of hospital catchment area definitions on hospital admissions forecasts: COVID-19 in England, September 2020-April 2021.

BACKGROUND: Defining healthcare facility catchment areas is a key step in predicting future healthcare demand in epidemic settings. Forecasts of hospitalisations can be informed by leading indicators measured at the community level. However, this relies on the definition of so-called catchment areas or the geographies whose populations make up the patients admitted to a given hospital, which are often not well-defined. Little work has been done to quantify the impact of hospital catchment area definitions on healthcare demand forecasting. METHODS: We made forecasts of local-level hospital admissions using a scaled convolution of local cases (as defined by the hospital catchment area) and delay distribution. Hospital catchment area definitions were derived from either simple heuristics (in which people are admitted to their nearest hospital or any nearby hospital) or historical admissions data (all emergency or elective admissions in 2019, or COVID-19 admissions), plus a marginal baseline definition based on the distribution of all hospital admissions. We evaluated predictive performance using each hospital catchment area definition using the weighted interval score and considered how this changed by the length of the predictive horizon, the date on which the forecast was made, and by location. We also considered the change, if any, in the relative performance of each definition in retrospective vs. real-time settings, or at different spatial scales. RESULTS: The choice of hospital catchment area definition affected the accuracy of hospital admission forecasts. The definition based on COVID-19 admissions data resulted in the most accurate forecasts at both a 7- and 14-day horizon and was one of the top two best-performing definitions across forecast dates and locations. The "nearby" heuristic also performed well, but less consistently than the COVID-19 data definition. The marginal distribution baseline, which did not include any spatial information, was the lowest-ranked definition. The relative performance of the definitions was larger when using case forecasts compared to future observed cases. All results were consistent across spatial scales of the catchment area definitions. CONCLUSIONS: Using catchment area definitions derived from context-specific data can improve local-level hospital admission forecasts. Where context-specific data is not available, using catchment areas defined by carefully chosen heuristics is a sufficiently good substitute. There is clear value in understanding what drives local admissions patterns, and further research is needed to understand the impact of different catchment area definitions on forecast performance where case trends are more heterogeneous.

Evaluating the use of social contact data to produce age-specific short-term forecasts of SARS-CoV-2 incidence in England.

Mathematical and statistical models can be used to make predictions of how epidemics may progress in the near future and form a central part of outbreak mitigation and control. Renewal equation based models allow inference of epidemiological parameters from historical data and forecast future epidemic dynamics without requiring complex mechanistic assumptions. However, these models typically ignore interaction between age groups, partly due to challenges in parameterising a time varying interaction matrix. Social contact data collected regularly during the COVID-19 epidemic provide a means to inform interaction between age groups in real-time. We developed an age-specific forecasting framework and applied it to two age-stratified time-series: incidence of SARS-CoV-2 infection, estimated from a national infection and antibody prevalence survey; and, reported cases according to the UK national COVID-19 dashboard. Jointly fitting our model to social contact data from the CoMix study, we inferred a time-varying next generation matrix which we used to project infections and cases in the four weeks following each of 29 forecast dates between October 2020 and November 2021. We evaluated the forecasts using proper scoring rules and compared performance with three other models with alternative data and specifications alongside two naive baseline models. Overall, incorporating age interaction improved forecasts of infections and the CoMix-data-informed model was the best performing model at time horizons between two and four weeks. However, this was not true when forecasting cases. We found that age group interaction was most important for predicting cases in children and older adults. The contact-data-informed models performed best during the winter months of 2020-2021, but performed comparatively poorly in other periods. We highlight challenges regarding the incorporation of contact data in forecasting and offer proposals as to how to extend and adapt our approach, which may lead to more successful forecasts in future.

Predictive performance of multi-model ensemble forecasts of COVID-19 across European nations.

Background: Short-term forecasts of infectious disease burden can contribute to situational awareness and aid capacity planning. Based on best practice in other fields and recent insights in infectious disease epidemiology, one can maximise the predictive performance of such forecasts if multiple models are combined into an ensemble. Here, we report on the performance of ensembles in predicting COVID-19 cases and deaths across Europe between 08 March 2021 and 07 March 2022. Methods: We used open-source tools to develop a public European COVID-19 Forecast Hub. We invited groups globally to contribute weekly forecasts for COVID-19 cases and deaths reported by a standardised source for 32 countries over the next 1-4 weeks. Teams submitted forecasts from March 2021 using standardised quantiles of the predictive distribution. Each week we created an ensemble forecast, where each predictive quantile was calculated as the equally-weighted average (initially the mean and then from 26th July the median) of all individual models' predictive quantiles. We measured the performance of each model using the relative Weighted Interval Score (WIS), comparing models' forecast accuracy relative to all other models. We retrospectively explored alternative methods for ensemble forecasts, including weighted averages based on models' past predictive performance. Results: Over 52 weeks, we collected forecasts from 48 unique models. We evaluated 29 models' forecast scores in comparison to the ensemble model. We found a weekly ensemble had a consistently strong performance across countries over time. Across all horizons and locations, the ensemble performed better on relative WIS than 83% of participating models' forecasts of incident cases (with a total N=886 predictions from 23 unique models), and 91% of participating models' forecasts of deaths (N=763 predictions from 20 models). Across a 1-4 week time horizon, ensemble performance declined with longer forecast periods when forecasting cases, but remained stable over 4 weeks for incident death forecasts. In every forecast across 32 countries, the ensemble outperformed most contributing models when forecasting either cases or deaths, frequently outperforming all of its individual component models. Among several choices of ensemble methods we found that the most influential and best choice was to use a median average of models instead of using the mean, regardless of methods of weighting component forecast models. Conclusions: Our results support the use of combining forecasts from individual models into an ensemble in order to improve predictive performance across epidemiological targets and populations during infectious disease epidemics. Our findings further suggest that median ensemble methods yield better predictive performance more than ones based on means. Our findings also highlight that forecast consumers should place more weight on incident death forecasts than incident case forecasts at forecast horizons greater than 2 weeks. Funding: AA, BH, BL, LWa, MMa, PP, SV funded by National Institutes of Health (NIH) Grant 1R01GM109718, NSF BIG DATA Grant IIS-1633028, NSF Grant No.: OAC-1916805, NSF Expeditions in Computing Grant CCF-1918656, CCF-1917819, NSF RAPID CNS-2028004, NSF RAPID OAC-2027541, US Centers for Disease Control and Prevention 75D30119C05935, a grant from Google, University of Virginia Strategic Investment Fund award number SIF160, Defense Threat Reduction Agency (DTRA) under Contract No. HDTRA1-19-D-0007, and respectively Virginia Dept of Health Grant VDH-21-501-0141, VDH-21-501-0143, VDH-21-501-0147, VDH-21-501-0145, VDH-21-501-0146, VDH-21-501-0142, VDH-21-501-0148. AF, AMa, GL funded by SMIGE - Modelli statistici inferenziali per governare l'epidemia, FISR 2020-Covid-19 I Fase, FISR2020IP-00156, Codice Progetto: PRJ-0695. AM, BK, FD, FR, JK, JN, JZ, KN, MG, MR, MS, RB funded by Ministry of Science and Higher Education of Poland with grant 28/WFSN/2021 to the University of Warsaw. BRe, CPe, JLAz funded by Ministerio de Sanidad/ISCIII. BT, PG funded by PERISCOPE European H2020 project, contract number 101016233. CP, DL, EA, MC, SA funded by European Commission - Directorate-General for Communications Networks, Content and Technology through the contract LC-01485746, and Ministerio de Ciencia, Innovacion y Universidades and FEDER, with the project PGC2018-095456-B-I00. DE., MGu funded by Spanish Ministry of Health / REACT-UE (FEDER). DO, GF, IMi, LC funded by Laboratory Directed Research and Development program of Los Alamos National Laboratory (LANL) under project number 20200700ER. DS, ELR, GG, NGR, NW, YW funded by National Institutes of General Medical Sciences (R35GM119582; the content is solely the responsibility of the authors and does not necessarily represent the official views of NIGMS or the National Institutes of Health). FB, FP funded by InPresa, Lombardy Region, Italy. HG, KS funded by European Centre for Disease Prevention and Control. IV funded by Agencia de Qualitat i Avaluacio Sanitaries de Catalunya (AQuAS) through contract 2021-021OE. JDe, SMo, VP funded by Netzwerk Universitatsmedizin (NUM) project egePan (01KX2021). JPB, SH, TH funded by Federal Ministry of Education and Research (BMBF; grant 05M18SIA). KH, MSc, YKh funded by Project SaxoCOV, funded by the German Free State of Saxony. Presentation of data, model results and simulations also funded by the NFDI4Health Task Force COVID-19 (https://www.nfdi4health.de/task-force-covid-19-2) within the framework of a DFG-project (LO-342/17-1). LP, VE funded by Mathematical and Statistical modelling project (MUNI/A/1615/2020), Online platform for real-time monitoring, analysis and management of epidemic situations (MUNI/11/02202001/2020); VE also supported by RECETOX research infrastructure (Ministry of Education, Youth and Sports of the Czech Republic: LM2018121), the CETOCOEN EXCELLENCE (CZ.02.1.01/0.0/0.0/17-043/0009632), RECETOX RI project (CZ.02.1.01/0.0/0.0/16-013/0001761). NIB funded by Health Protection Research Unit (grant code NIHR200908). SAb, SF funded by Wellcome Trust (210758/Z/18/Z).

The association of COVID-19 incidence with temperature, humidity, and UV radiation - A global multi-city analysis.

BACKGROUND AND AIM: The associations between COVID-19 transmission and meteorological factors are scientifically debated. Several studies have been conducted worldwide, with inconsistent findings. However, often these studies had methodological issues, e.g., did not exclude important confounding factors, or had limited geographic or temporal resolution. Our aim was to quantify associations between temporal variations in COVID-19 incidence and meteorological variables globally. METHODS: We analysed data from 455 cities across 20 countries from 3 February to 31 October 2020. We used a time-series analysis that assumes a quasi-Poisson distribution of the cases and incorporates distributed lag non-linear modelling for the exposure associations at the city-level while considering effects of autocorrelation, long-term trends, and day of the week. The confounding by governmental measures was accounted for by incorporating the Oxford Governmental Stringency Index. The effects of daily mean air temperature, relative and absolute humidity, and UV radiation were estimated by applying a meta-regression of local estimates with multi-level random effects for location, country, and climatic zone. RESULTS: We found that air temperature and absolute humidity influenced the spread of COVID-19 over a lag period of 15 days. Pooling the estimates globally showed that overall low temperatures (7.5 °C compared to 17.0 °C) and low absolute humidity (6.0 g/m3 compared to 11.0 g/m3) were associated with higher COVID-19 incidence (RR temp =1.33 with 95%CI: 1.08; 1.64 and RR AH =1.33 with 95%CI: 1.12; 1.57). RH revealed no significant trend and for UV some evidence of a positive association was found. These results were robust to sensitivity analysis. However, the study results also emphasise the heterogeneity of these associations in different countries. CONCLUSION: Globally, our results suggest that comparatively low temperatures and low absolute humidity were associated with increased risks of COVID-19 incidence. However, this study underlines regional heterogeneity of weather-related effects on COVID-19 transmission.

Effect of Fatigue on Health-Related Quality of Life and Work Productivity in Psoriatic Arthritis: Findings From a Real-World Survey.

OBJECTIVE: To evaluate fatigue frequency and severity among patients with psoriatic arthritis (PsA) and assess the effect of fatigue severity on patient-reported outcome measures (PROMs) assessing quality of life, function, and work productivity. METHODS: Data were derived from the Adelphi Disease Specific Programme, a cross-sectional survey conducted in 2018 in the United States and Europe. Patients had physician-confirmed PsA. Fatigue was collected as a binary variable and through its severity (0-10 scale, using the 12-item Psoriatic Arthritis Impact of Disease fatigue question) from patients; physicians also reported patient fatigue (yes/no). Other PROMs included the 5-level EuroQol 5-dimension questionnaire (EQ-5D-5L) for health-related quality of life (HRQOL), Health Assessment Questionnaire-Disability Index, and Work Productivity and Activity Impairment Questionnaire. Multivariate linear regression was used to evaluate the association between fatigue severity and other PROMs. RESULTS: Among the 831 included patients (mean age 47.5 yrs, mean disease duration 5.3 yrs, 46.9% female, 48.1% receiving a biologic), fatigue was reported by 78.3% of patients. Patients with greater fatigue severity had greater disease duration, PsA severity, pain levels, body surface area affected by psoriasis, and swollen and tender joint counts (all P < 0.05). In multivariate analyses, patients with greater fatigue severity experienced worse physical functioning, HRQOL, and work productivity (all P < 0.001). Presence of fatigue was underreported by physicians (reported in only 32% of patients who self-reported fatigue). CONCLUSION: Prevalence of patient-reported fatigue was high among patients with PsA and underrecognized by physicians. Fatigue severity was associated with altered physical functioning, work productivity, and HRQOL.

Patient-psychiatrist discordance and drivers of prescribing long-acting injectable antipsychotics for schizophrenia management in the real-world: a point-in-time survey.

BACKGROUND: To assess discordance between psychiatrists and their patients with schizophrenia regarding disease management and understand drivers of prescribing long-acting injectable (LAI) antipsychotics. METHODS: Data were collected via the Adelphi Schizophrenia Disease Specific Programme™, a point-in-time real-world international survey of psychiatrists and their consulting patients with schizophrenia, conducted in 2019. Psychiatrists completed an attitudinal survey on schizophrenia management and provided patient profiles for their next 10 adult consulting patients. The same patients voluntarily completed patient self-completion forms. Disease severity and improvement were assessed via physician-reported Clinical Global Impression scale; patients' adherence to treatment was rated through a 3-point scale (1=not at all adherent, 3=fully adherent). RESULTS: Four hundred sixty-six psychiatrists provided data for 4345 patients (1132 receiving a LAI; 3105 on non-LAI treatment; 108 not on treatment). LAIs were more commonly prescribed to patients with severe schizophrenia, with varying reasons for prescribing. Globally, only slight agreement was observed between psychiatrists and patients for Clinical Global Impression severity of illness (κ=0.174) and level of improvement on treatment (κ=0.204). There was moderate agreement regarding level of adherence to treatment (κ=0.524). Reasons for non-adherence did not reach a level of agreement greater than fair. CONCLUSIONS: Our real-world survey found that LAIs were more often reserved for severe schizophrenia patients and improving adherence was a key driver for prescribing. However, compared with the patients themselves, psychiatrists tended to underestimate patients' disease severity and overestimate their adherence.

Comparative assessment of methods for short-term forecasts of COVID-19 hospital admissions in England at the local level.

BACKGROUND: Forecasting healthcare demand is essential in epidemic settings, both to inform situational awareness and facilitate resource planning. Ideally, forecasts should be robust across time and locations. During the COVID-19 pandemic in England, it is an ongoing concern that demand for hospital care for COVID-19 patients in England will exceed available resources. METHODS: We made weekly forecasts of daily COVID-19 hospital admissions for National Health Service (NHS) Trusts in England between August 2020 and April 2021 using three disease-agnostic forecasting models: a mean ensemble of autoregressive time series models, a linear regression model with 7-day-lagged local cases as a predictor, and a scaled convolution of local cases and a delay distribution. We compared their point and probabilistic accuracy to a mean-ensemble of them all and to a simple baseline model of no change from the last day of admissions. We measured predictive performance using the weighted interval score (WIS) and considered how this changed in different scenarios (the length of the predictive horizon, the date on which the forecast was made, and by location), as well as how much admissions forecasts improved when future cases were known. RESULTS: All models outperformed the baseline in the majority of scenarios. Forecasting accuracy varied by forecast date and location, depending on the trajectory of the outbreak, and all individual models had instances where they were the top- or bottom-ranked model. Forecasts produced by the mean-ensemble were both the most accurate and most consistently accurate forecasts amongst all the models considered. Forecasting accuracy was improved when using future observed, rather than forecast, cases, especially at longer forecast horizons. CONCLUSIONS: Assuming no change in current admissions is rarely better than including at least a trend. Using confirmed COVID-19 cases as a predictor can improve admissions forecasts in some scenarios, but this is variable and depends on the ability to make consistently good case forecasts. However, ensemble forecasts can make forecasts that make consistently more accurate forecasts across time and locations. Given minimal requirements on data and computation, our admissions forecasting ensemble could be used to anticipate healthcare needs in future epidemic or pandemic settings.

Physician-reported Clinical Unmet Needs, Burden and Treatment Patterns of Paediatric Psoriasis Patients: A US and EU Real-world Evidence Study.

This study is a retrospective analysis using data collected from the Adelphi Paediatric Psoriasis Disease-Specific Programme cross-sectional survey. Despite being treated for their psoriasis, a substantial proportion of paediatric patients presented with moderate (18.3%) or severe (1.3%) disease at sampling; 42.9% and 92.0% had a body surface area (BSA) of >10%, and 38.8% and 100.0% had a Psoriasis Area Severity Index (PASI) score >10, respectively. Overall, 69.9% of patients had only ever been treated with a topical therapy for their psoriasis. For patients with moderate or severe disease at sampling, 16.3% and 14.4% were currently receiving conventional systemics or biologic therapy, respectively. There is a clinical unmet need in this paediatric population; a considerable percentage of patients still experienced moderate or severe disease and persistent psoriasis symptoms, with numerous body areas affected. A significant proportion of patients were undertreated, which may explain the high burden of disease observed.

Comparative assessment of methods for short-term forecasts of COVID-19 hospital admissions in England at the local level.

BACKGROUND: Forecasting healthcare demand is essential in epidemic settings, both to inform situational awareness and facilitate resource planning. Ideally, forecasts should be robust across time and locations. During the COVID-19 pandemic in England, it is an ongoing concern that demand for hospital care for COVID-19 patients in England will exceed available resources. METHODS: We made weekly forecasts of daily COVID-19 hospital admissions for National Health Service (NHS) Trusts in England between August 2020 and April 2021 using three disease-agnostic forecasting models: a mean ensemble of autoregressive time series models, a linear regression model with 7-day-lagged local cases as a predictor, and a scaled convolution of local cases and a delay distribution. We compared their point and probabilistic accuracy to a mean-ensemble of them all, and to a simple baseline model of no change from the last day of admissions. We measured predictive performance using the Weighted Interval Score (WIS) and considered how this changed in different scenarios (the length of the predictive horizon, the date on which the forecast was made, and by location), as well as how much admissions forecasts improved when future cases were known. RESULTS: All models outperformed the baseline in the majority of scenarios. Forecasting accuracy varied by forecast date and location, depending on the trajectory of the outbreak, and all individual models had instances where they were the top- or bottom-ranked model. Forecasts produced by the mean-ensemble were both the most accurate and most consistently accurate forecasts amongst all the models considered. Forecasting accuracy was improved when using future observed, rather than forecast, cases, especially at longer forecast horizons. CONCLUSIONS: Assuming no change in current admissions is rarely better than including at least a trend. Using confirmed COVID-19 cases as a predictor can improve admissions forecasts in some scenarios, but this is variable and depends on the ability to make consistently good case forecasts. However, ensemble forecasts can make forecasts that make consistently more accurate forecasts across time and locations. Given minimal requirements on data and computation, our admissions forecasting ensemble could be used to anticipate healthcare needs in future epidemic or pandemic settings.

Correction: Practical considerations for measuring the effective reproductive number, Rt.

[This corrects the article DOI: 10.1371/journal.pcbi.1008409.].

Associations Between Relapses and Psychosocial Outcomes in Patients With Schizophrenia in Real-World Settings in the United States.

Aim: To assess associations between relapses and psychosocial outcomes in adult patients with schizophrenia treated in United States (US) healthcare settings. Methods: Data were derived from a point-in-time survey of psychiatrists and their patients with schizophrenia conducted across the US, France, Spain, China, and Japan between July and October 2019. For the purposes of this analysis, only data from US practitioners and patients were included. Disease-specific programmes (DSPs) are large surveys with a validated methodology conducted in clinical practise; they describe current disease management, disease burden, and associated treatment effects (clinical and physician-perceived). Participating psychiatrists completed patient record forms for their next 10 consecutive adult consulting patients with schizophrenia, with the same patients invited to voluntarily complete a patient self-completion (PSC) questionnaire. Surveys contained questions on the patients' disease background, treatment history, prior hospitalisation due to schizophrenia relapse and a series of psychosocial outcomes. Associations between relapses in the last 12 months and psychosocial outcomes were examined using multiple regression. Results: A total of 124 psychiatrists provided data on 1,204 patients. Of these, 469 patients (mean age, 39.6 years; 56.5% male) had known hospitalisation history for the last 12 months and completed a PSC; 116 (24.7%) patients had ≥1 relapse. Compared to patients without relapses, patients who relapsed were more likely to be homeless, unemployed, previously incarcerated, and currently have difficulties living independently (all p < 0.05). Patients who experience a relapse also had greater working impairment and poorer quality of life compared with those who did not relapse. In general, psychosocial outcomes became poorer with an increasing number of relapses. Conclusions: In this population of patients with schizophrenia from the US, relapse was significantly associated with poor psychosocial outcomes, with a greater number of relapses predicting worse outcomes. Early intervention to reduce the risk of relapse may improve psychosocial outcomes in patients with schizophrenia.

A cross-sectional analysis of meteorological factors and SARS-CoV-2 transmission in 409 cities across 26 countries.

There is conflicting evidence on the influence of weather on COVID-19 transmission. Our aim is to estimate weather-dependent signatures in the early phase of the pandemic, while controlling for socio-economic factors and non-pharmaceutical interventions. We identify a modest non-linear association between mean temperature and the effective reproduction number (Re) in 409 cities in 26 countries, with a decrease of 0.087 (95% CI: 0.025; 0.148) for a 10 °C increase. Early interventions have a greater effect on Re with a decrease of 0.285 (95% CI 0.223; 0.347) for a 5th - 95th percentile increase in the government response index. The variation in the effective reproduction number explained by government interventions is 6 times greater than for mean temperature. We find little evidence of meteorological conditions having influenced the early stages of local epidemics and conclude that population behaviour and government interventions are more important drivers of transmission.

How to improve outbreak response: a case study of integrated outbreak analytics from Ebola in Eastern Democratic Republic of the Congo.

The emerging field of outbreak analytics calls attention to the need for data from multiple sources to inform evidence-based decision making in managing infectious diseases outbreaks. To date, these approaches have not systematically integrated evidence from social and behavioural sciences. During the 2018-2020 Ebola outbreak in Eastern Democratic Republic of the Congo, an innovative solution to systematic and timely generation of integrated and actionable social science evidence emerged in the form of the Cellulle d'Analyse en Sciences Sociales (Social Sciences Analytics Cell) (CASS), a social science analytical cell. CASS worked closely with data scientists and epidemiologists operating under the Epidemiological Cell to produce integrated outbreak analytics (IOA), where quantitative epidemiological analyses were complemented by behavioural field studies and social science analyses to help better explain and understand drivers and barriers to outbreak dynamics. The primary activity of the CASS was to conduct operational social science analyses that were useful to decision makers. This included ensuring that research questions were relevant, driven by epidemiological data from the field, that research could be conducted rapidly (ie, often within days), that findings were regularly and systematically presented to partners and that recommendations were co-developed with response actors. The implementation of the recommendations based on CASS analytics was also monitored over time, to measure their impact on response operations. This practice paper presents the CASS logic model, developed through a field-based externally led consultation, and documents key factors contributing to the usefulness and adaption of CASS and IOA to guide replication for future outbreaks.

Exploring surveillance data biases when estimating the reproduction number: with insights into subpopulation transmission of COVID-19 in England.

The time-varying reproduction number (Rt: the average number of secondary infections caused by each infected person) may be used to assess changes in transmission potential during an epidemic. While new infections are not usually observed directly, they can be estimated from data. However, data may be delayed and potentially biased. We investigated the sensitivity of Rt estimates to different data sources representing COVID-19 in England, and we explored how this sensitivity could track epidemic dynamics in population sub-groups. We sourced public data on test-positive cases, hospital admissions and deaths with confirmed COVID-19 in seven regions of England over March through August 2020. We estimated Rt using a model that mapped unobserved infections to each data source. We then compared differences in Rt with the demographic and social context of surveillance data over time. Our estimates of transmission potential varied for each data source, with the relative inconsistency of estimates varying across regions and over time. Rt estimates based on hospital admissions and deaths were more spatio-temporally synchronous than when compared to estimates from all test positives. We found these differences may be linked to biased representations of subpopulations in each data source. These included spatially clustered testing, and where outbreaks in hospitals, care homes, and young age groups reflected the link between age and severity of the disease. We highlight that policy makers could better target interventions by considering the source populations of Rt estimates. Further work should clarify the best way to combine and interpret Rt estimates from different data sources based on the desired use. This article is part of the theme issue 'Modelling that shaped the early COVID-19 pandemic response in the UK'.

Real-time monitoring of COVID-19 dynamics using automated trend fitting and anomaly detection.

As several countries gradually release social distancing measures, rapid detection of new localized COVID-19 hotspots and subsequent intervention will be key to avoiding large-scale resurgence of transmission. We introduce ASMODEE (automatic selection of models and outlier detection for epidemics), a new tool for detecting sudden changes in COVID-19 incidence. Our approach relies on automatically selecting the best (fitting or predicting) model from a range of user-defined time series models, excluding the most recent data points, to characterize the main trend in an incidence. We then derive prediction intervals and classify data points outside this interval as outliers, which provides an objective criterion for identifying departures from previous trends. We also provide a method for selecting the optimal breakpoints, used to define how many recent data points are to be excluded from the trend fitting procedure. The analysis of simulated COVID-19 outbreaks suggests ASMODEE compares favourably with a state-of-art outbreak-detection algorithm while being simpler and more flexible. As such, our method could be of wider use for infectious disease surveillance. We illustrate ASMODEE using publicly available data of National Health Service (NHS) Pathways reporting potential COVID-19 cases in England at a fine spatial scale, showing that the method would have enabled the early detection of the flare-ups in Leicester and Blackburn with Darwen, two to three weeks before their respective lockdown. ASMODEE is implemented in the free R package trendbreaker. This article is part of the theme issue 'Modelling that shaped the early COVID-19 pandemic response in the UK'.

Implications of the school-household network structure on SARS-CoV-2 transmission under school reopening strategies in England.

In early 2020 many countries closed schools to mitigate the spread of SARS-CoV-2. Since then, governments have sought to relax the closures, engendering a need to understand associated risks. Using address records, we construct a network of schools in England connected through pupils who share households. We evaluate the risk of transmission between schools under different reopening scenarios. We show that whilst reopening select year-groups causes low risk of large-scale transmission, reopening secondary schools could result in outbreaks affecting up to 2.5 million households if unmitigated, highlighting the importance of careful monitoring and within-school infection control to avoid further school closures or other restrictions.

Practical considerations for measuring the effective reproductive number, Rt.

Estimation of the effective reproductive number Rt is important for detecting changes in disease transmission over time. During the Coronavirus Disease 2019 (COVID-19) pandemic, policy makers and public health officials are using Rt to assess the effectiveness of interventions and to inform policy. However, estimation of Rt from available data presents several challenges, with critical implications for the interpretation of the course of the pandemic. The purpose of this document is to summarize these challenges, illustrate them with examples from synthetic data, and, where possible, make recommendations. For near real-time estimation of Rt, we recommend the approach of Cori and colleagues, which uses data from before time t and empirical estimates of the distribution of time between infections. Methods that require data from after time t, such as Wallinga and Teunis, are conceptually and methodologically less suited for near real-time estimation, but may be appropriate for retrospective analyses of how individuals infected at different time points contributed to the spread. We advise caution when using methods derived from the approach of Bettencourt and Ribeiro, as the resulting Rt estimates may be biased if the underlying structural assumptions are not met. Two key challenges common to all approaches are accurate specification of the generation interval and reconstruction of the time series of new infections from observations occurring long after the moment of transmission. Naive approaches for dealing with observation delays, such as subtracting delays sampled from a distribution, can introduce bias. We provide suggestions for how to mitigate this and other technical challenges and highlight open problems in Rt estimation.

Practical considerations for measuring the effective reproductive number, R t.

Estimation of the effective reproductive number, R t , is important for detecting changes in disease transmission over time. During the COVID-19 pandemic, policymakers and public health officials are using R t to assess the effectiveness of interventions and to inform policy. However, estimation of R t from available data presents several challenges, with critical implications for the interpretation of the course of the pandemic. The purpose of this document is to summarize these challenges, illustrate them with examples from synthetic data, and, where possible, make recommendations. For near real-time estimation of R t , we recommend the approach of Cori et al. (2013), which uses data from before time t and empirical estimates of the distribution of time between infections. Methods that require data from after time t, such as Wallinga and Teunis (2004), are conceptually and methodologically less suited for near real-time estimation, but may be appropriate for retrospective analyses of how individuals infected at different time points contributed to spread. We advise against using methods derived from Bettencourt and Ribeiro (2008), as the resulting R t estimates may be biased if the underlying structural assumptions are not met. Two key challenges common to all approaches are accurate specification of the generation interval and reconstruction of the time series of new infections from observations occurring long after the moment of transmission. Naive approaches for dealing with observation delays, such as subtracting delays sampled from a distribution, can introduce bias. We provide suggestions for how to mitigate this and other technical challenges and highlight open problems in R t estimation.

Correlations between stochastic endemic infection in multiple interacting subpopulations.

Heterogeneity plays an important role in the emergence, persistence and control of infectious diseases. Metapopulation models are often used to describe spatial heterogeneity, and the transition from random- to heterogeneous-mixing is made by incorporating the interaction, or coupling, within and between subpopulations. However, such couplings are difficult to measure explicitly; instead, their action through the correlations between subpopulations is often all that can be observed. We use moment-closure methods to investigate how the coupling and resulting correlation are related, considering systems of multiple identical interacting populations on highly symmetric complex networks: the complete network, the k-regular tree network, and the star network. We show that the correlation between the prevalence of infection takes a relatively simple form and can be written in terms of the coupling, network parameters and epidemiological parameters only. These results provide insight into the effect of metapopulation network structure on endemic disease dynamics, and suggest that detailed case-reporting data alone may be sufficient to infer the strength of between population interaction and hence lead to more accurate mathematical descriptions of infectious disease behaviour.

Correlations between stochastic epidemics in two interacting populations.

It is increasingly apparent that heterogeneity in the interaction between individuals plays an important role in the dynamics, persistence, evolution and control of infectious diseases. In epidemic modelling two main forms of heterogeneity are commonly considered: spatial heterogeneity due to the segregation of populations and heterogeneity in risk at the same location. The transition from random-mixing to heterogeneous-mixing models is made by incorporating the interaction, or coupling, within and between subpopulations. However, such couplings are difficult to measure explicitly; instead, their action through the correlations between subpopulations is often all that can be observed. Here, using moment-closure methodology supported by stochastic simulation, we investigate how the coupling and resulting correlation are related. We focus on the simplest case of interactions, two identical coupled populations, and show that for a wide range of parameters the correlation between the prevalence of infection takes a relatively simple form. In particular, the correlation can be approximated by a logistic function of the between population coupling, with the free parameter determined analytically from the epidemiological parameters. These results suggest that detailed case-reporting data alone may be sufficient to infer the strength of between population interaction and hence lead to more accurate mathematical descriptions of infectious disease behaviour.